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Abstract 

Measurements  of  heat  transfer  parameters  in  pulsating  and  steady, 
forced,  turbulent  air  flow  through  a  vertical  circular  tube,  which  was 
resistively  heated,  are  reported.  Inlet  Reynolds  numbers  varied  from 
19,000  to  102,000;  Mach  numbers  were  0.15  or  below;  pulsation  frequencies 
ranged  from  2.1  to  3.6  Hz  and  the  peak-to-peak  pressure  fluctuations 
varied  from  9  to  29  percent  of  the  mean  pressure.  At  these  conditions 
the  non-dimensional  frequency  a  *  r/2irf/v,  varied  from  about  4  to  7  1/2; 
in  laminar  flow,  quasi-steady  approximations  became  weak  when  this 
frequency  becomes  greater  than  about  two,  but  for  turbulent  heat  trans¬ 
fer  in  this  Reynolds  number  range  the  limitations  are  still  to  be  deter¬ 
mined.  Direct  comparison  between  steady  and  pulsating  runs  at  the  same 
values  of  the  control  parameters  essentially  confirmed  the  quasi-steady 
analyses  for  Re  x  5  x  10  . 
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NOMENCLATURE 


Definition _ 

Speed  of  sound 

Specific  heat  at  constant  pressure 

Diameter  of  tube 

Frequency,  cycles/ sec. 

Average  mass  flux,  pV 

4  2 

Grashof  number,  gD  q"/ (\<  kT) 

w 

Gravitational  acceleration 

2 

Dimensional  conversion  factor,  e.g.,  32.174  lbm-ft/ (lbf ‘sec  ) 

Heat  transfer  coefficient 

Thermal  conductivity 

Length 

Mach  number 

Exponent  in  heat  transfer  correlation 

Order  of 

Pressure 

Mean  pressure 

Period 

Prandtl  number,  yc^/p 
Non-dimensional  heat  flux,  q^V (pVc^T^) 

Radius  of  tube 

Reynolds  number,  pVL/y  (or,  for  tube,  4,,/'n’Dy 

m 

Strouhal  number 
Temperature 
Bulk  velocity 


NOMENCLATURE—Continued 


Definition _ 

Mean  bulk  velocity 
Non-dimensional  velocity  amplitude 

Subscripts 

Bulk 

Inlet 

Mean 

Pulsed 

Steady 

Wall 

Greek  symbols 

Non-dimensional  frequency,  r/w/V 
Time 

Wave  length 
Dynamic  viscoslL/ 

Kinematic  viscosity,  y/p 
Density 

Angular  frequency,  2irf 


1.  INTRODUCTION 


Current  Naval  propulsion  plants  are  powered  by  variations  of  the 
Rankine  cycle  (steam)  or  the  open  gas  turbine  cycle  (air  and  combustion 
products)  plus  some  diesel  engines  in  small  ships.  Alternative  power 
systems  suggested  include  the  closed  gas  turbine  cycle  and  cycles  in¬ 
volving  dissociation  of  the  working  fluid  in  either  a  Rankine  or  a  gas 
cycle.  These  latter  two  are  believed  to  offer  the  potential  of  sub¬ 
stantial  improvement  in  the  power-to-weight  ratio  of  the  propulsion 
plant. 

Convective  heat  transfer  provides  the  dominant  thermal  resistance 
in  several  components  of  conventional  steam  power  plants  as  well  as  in 
all  heat  transfer  components  in  gaseous  cycles.  In  order  to  design 
propulsion  plants  with  better  power-to-weight  ratios  and/or  fuel  econ¬ 
omy  reliably,  the  convective  heat  transfer  must  be  predicted  accurately. 
The  present  studies  consider  basic  problems  in  convective  heat  transfer 
and  flow  friction  that  are  important  in  all  of  the  above. 

Pulsating  Internal  gas  flows  occur  frequently  in  practise.  In 
some  cases  the  pulsations  are  unintentional  consequences  of  the  flow 
equipment  employed,  e.g. ,  reciprocating  compressors,  and  in  others  they 
are  meant  to  accomplish  a  purpose  such  as  improving  convective  heat  trans¬ 
fer  [Lemlich,  1961].  In  either  case  it  is  Important  for  the  designer  to 
know  the  effects  of  the  pulsations  on  the  time-mean  heat  transfer  param¬ 
eters  for  his/her  proposed  range  of  operations.  Also,  the  so-called  mass 


2 


flowmeter,  or  thermal  flowmeter,  employs  wall  temperature  measurement 
for  convective  heat  transfer  in  a  tube  in  order  to  deduce  the  mass  flow 
rate  of  a  gas  [Zinsmeister  and  Dixon,  1966];  the  extent  to  which  its 
steady  flow  calibration  is  valid  for  non-steady  flows  also  depends  on 

the  magnitude  of  the  effects  of  pulsating  flows. 

The  present  study  utilizes  a  vertical  circular  tube  with  flow  pulsa¬ 
tions  superposed  on  the  throughflow.  The  frequencies  are  sufficiently 
low  that  possible  acoustical  resonances  are  not  expected  to  be  important. 


1.1.  Governing  parameters  of  non-steady  flow 

Smolderen  [1977]  provides  a  general  introduction  to  the  theory  of 
unsteady  fluid  dynamic  phenomena  and  identifies  pertinent  non-dimensional 
parameters.  The  Strouhal  number 


Str. 


JL 

VP 


fL 

V 


is  the  ratio  of  the  length  of  the  flow  path  to  the  distance  traveled  by  a 
typical  fluid  particle  during  the  period  of  a  cycle.  It  is  also  called 
the  reduced  frequency  by  some  authors.  If  Str^  «  1,  unsteady  terms  in 
the  governing  equations  become  negligible  and  the  flow  can  be  treated 
as  quasi-steady.  For  acoustic  disturbances,  the  ratio  of  the  length  of 
the  flow  path  to  the  wave  length  of  sound  waves, 


X  a  P 
o 


^  ■  H*Str 

A 


L 


o 


is  an  important  parameter.  If  it  is  small,  the  flow  is  considered  quasi¬ 
steady. 
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For  consideration  of  viscous  effects,  one  examines  the  ratio 
of  the  length  to  the  penetration  depth  for  vorticity  diffusion,  which  is 
0 (/jp)  in  viscous  flows.  The  following  equivalency  results: 


RVSttL 


If  this  grouping  is  small,  vorticity  can  diffuse  through  the  entire 
field  in  a  fraction  of  a  period.  For  fully  developed  flow  in  a  tube, 
the  radius  is  the  logical  characteristic  length,  so  r/f/v  is  the  appro¬ 
priate  parameter;  many  authors  use  the  angular  frequency,  u  ■  2irf, 

and  for  this  non-dimensional  frequency  as  o  ■  r/co/v  [Baird  et  al.,  1966]. 
For  a  thermal  problem,  thermal  diffusion  may  becosie  important  and 

the  appropriate  length  scale  is  the  penetration  depth  of  a  thermal  dis¬ 
turbance,  0(/l kP/pCp).  The  related  parameter  becomes 


L2  vpcn  L2 

k*7pcp  k  VP  "  Pr,ReL,strL 


or,  for  fully  developed  flow  in  a  tube, 


kP/pc 


Pr 

2tt 


2 

a 


If  this  parameter  is  small,  the  thermal  disturbance  penetrates  the  entire 
field  and  the  thermal  treatment  can  be  quasi-steady  provided  the  velocity 
field  is  also  quasi-steady. 
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Richardson  [1967]  presents  additional  background  and  develops  comparable 
non-dimensional  parameters  via  a  different  approach.  For  turbulent  flows, 
he  suggests  that  an  "a.c.  boundary  layer  thickness,"  A>/u  (a  vorticity 
penetration  depth  from  Smolderen)  be  compared  to  the  "laminar  sublayer 
thickness,"  y+  %  5,  giving  5^uu7TJ”  as  a  parameter  and  a  criterion  for 
interation  with  non-steady  effects  in  turbulent  flows. 

1.2.  Related  work 

Lemlich  [1961],  Richardson  [1967],  Barnett  [1970]  and  others  present 
reviews  of  previous  studies  of  heat  transfer  to  pulsating  flows;  the  litera¬ 
ture  survey  which  accompanied  the  present  work  is  included  herein  as  Appendix  A. 

For  fully  developed,  laminar  flow  in  a  tube  Mull in  and  Created  [1980] 
point  out  that  for  low  values  of  the  viscous  frequency  parameter, 
a  -  r/U7v~ ,  quasi-steady  assumptions  should  hold,  while  at  a  >  10  the 
velocity  profile  changes  substantially. 

Lemlich  [1961]  shows  that  for  fully  developed  flow  with  constant 
properties,  if  the  heat  transfer  coefficient  varies  as  V°  for  steady 
flow,  the  quasi-steady  prediction  of  "enhancement"  would  be 


provided  the  flow  does  not  reverse.  He  notes  that  the  result  is  inde¬ 
pendent  of  frequency  and,  if  n  <  1,  a  decrease  in  h  will  be  predicted. 

Mueller  [1957]  calculated  the  enhancement  for  a  rectangular  wave 
representation  of  V(9)  and  showed  that  the  predicted  reduction  for 


dimensionless  amplitudes  of  1  and  less  (i.e.,  non-reversing)  would  be 
small.  For  a  symmetric  wave  V(6)  of  dimensionless  amplitude  0.5,  the 
reduction  would  be  about  2  to  3  percent.  He  also  derived  a  criterion 
for  quasi-steady  turbulent  flow  by  considering  the  "laminar"  sublayer 
and  presented  it  as  a  function  of  dimensionless  amplitude  and  the 
reciprocal  of  a^. 

Baird  et  al.  [1966]  predict  that  for  pulsations  of  the  form 

V  -  V  (1  +  Y  cos  u0 ) 
in 


the  enhancement  would  be 


<1  +  Y  cos  u>e) 


n 


de 


under  the  quasi-steady  assumptions.  This  value  is  negative  but  yields 
only  a  small  reduction  for  Y  <  1.  They  extend  the  treatment  to  reversing 
flows  and  show  that  for  Y  >  1.5  the  heat  transfer  coefficient  is  predicted 
to  improve  substantially.  For  their  conditions  they  estimate  the  limit 
of  the  quasi-steady  approximation  to  be  a  <  7. 4, based  on  residence  time 
considerations . 

Barnett  and  Vachon  [1970]  solved  the  non-steady  problem  approximately 
for  fully  developed, turbulent  pipe  flow.  The  turbulent  diffusivities 
were  taken  as  stationary  at  the  values  corresponding  to  the  mean  flow. 

They  predicted  significant  increases  of  heat  transfer  at  low  frequencies 
and  large  amplitudes;  the  calculated  Improvement  also  increases  as  the 
Prandtl  number  is  reduced. 


6 


Thomas  [1974]  applied  a  surface  renewal  model  and  predicted  only  slight 

changes  as  the  pulsation  frequency  increases  in  a  water  flow.  He  quotes 

2  ^  _ 

Lu  [1972]  and  Brown,  Margolis  and  Shah  [1969]  as  recommending  4a  •  0.1  Re^ 
as  a  limit  for  quasi-steady  flow. 

Measurements  of  heat  transfer  to  pulsating  air  flow  have  emphasized 
high  frequency  conditions  where  acoustical  resonances  become  Important 
[Lemlich,  1961].  Of  those  which  concern  low  frequency  pulsations  about 
half  the  experiments  used  steam-to-alr  heat  exchangers  which  do  not 
normally  permit  precise  comparisons.  Enhancements  and  reductions  of  the 
order  of  10  to  80  percent  have  been  reported  [Havemann  et  al. ,  1956; 
Chalitbhan,  1959]. 

With  electrical  heating  giving  an  approximately  constant  wall  heat 
flux,  Romle  [1956]  conducted  studies  at  Re  ■  5000  and  found  reductions  in 
heat  transfer  coefficients  of  up  to  10  percent  and  enhancements  up  to 
about  20  percent.  For  a  range  540  <  Re  <  11,000  in  comparable  apparatus, 
Mamayev,  Nosov  and  Syromyatnlkov  [1976]  found  apparent  enhancement  for 
low  frequencies  and  reductions  at  their  highest  frequency.  Data  at 
higher  Reynolds  numbers  do  not  appear  to  be  available  with  a  constant 
wall  heat  flux  as  the  boundary  condition  and  with  moderate  or  low  pulsation 
frequencies. 

In  summary,  the  analyses  generally  predict  only  a  slight  modification 
of  heat  transfer  parameters  in  pulsating  turbulent  flow  whereas  experi¬ 
ments  have  found  larger  effects.  Therefore,  it  is  important  to  measure  the 
heat  transfer  parameters  in  typical  flows  where  data  are  not  available 
in  order  to  test  the  analytical  predictions  for  normal  operating  con¬ 
ditions. 
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1.3.  Objective 

The  purpose  o£  the  present  study  Is  to  determine  by  measurement 
whether  the  predictions  of  quasi-steady  analyses  -  e.g.„  no  significant 
modification  of  heat  transfer  parameters  -  are  reasonable  for  pulsed 
turbulent  flow  at  typical  Reynolds  numbers  and  moderate  frequencies  and 
amplitudes.  In  order  to  avoid  some  uncertainties  of  previous  studies, 
a  direct  comparison  method  is  employed. 

2.  EXPERIMENT 
2.1.  Apparatus 

Data  were  obtained  with  upflow  through  a  vertical  circular  tube 
heated  reslstively.  The  apparatus  and  procedures  were  similar  to 
those  used  by  Pickett,  Taylor  and  McEligot  [1979]. 

Measurements  were  obtained  with  the  open  loop  apparatus  shown 
schematically  in  Figure  1.  A  regulated  gas  supply  flowed  directly 
into  the  system  or  through  a  single-acting  "Gas  Booster  Pump"  from 

Haskell  Manufacturing  Company.  The  first  source  was  used  for  steady 
flow  and  the  second  provided  pulsations  superposed  on  a  main  flow. 

From  the  pump  the  gas  flowed  through  3.2  meters  of  9.5mm  diameter 
tubing  to  a  heat  exchanger  of  14mm  ID  and  1.9m  long.  The  flow 
continued  through  the  same  diameter  tube  0.67m  to  a  fitting  con¬ 
taining  a  thermocouple  used  to  measure  the  gas  temperature  and  another 
fitting  in  which  a  Model  SCD  147  pressure  transducer  from  Data  Instru¬ 
ments,  Inc.  was  located  to  measure  pressure  fluctuations.  The  flow 
continued  through  the  same  diameter  tubing  for  another  7.6cm  and  at 


Figure  1.  Schematic  diagram  of  experimental  apparatus. 


that  point  the  flow  entered  a  tube  1.54m  long  with  a  5.87nsn  diameter; 
this  tube  served  as  an  alternate  flowmeter  for  steady  flow  runs.  At 
the  exit  of  this  tubular  flowmeter  was  an  elbow  which  served  as  the 
entrance  to  the  test  section. 

The  test  section  was  a  tube  of  Inconel  600  with  inside  diameter 
of  5.87mm  and  wall  thickness  of  0.292  mm.  It  consisted  of  a  heated 
section  of  60  diameters  in  length  preceded  by  an  unheated  section  of 
60  diameters  and  followed  by  an  unheated  section  of  56  diameters. 

A  Kulite  model  XT-140-100G  subminiature  pressure  transducer  was  mounted 
flush  witl.  the  inside  of  the  tube  immediately  beyond  the  elbow  at  the 
entry  of  the  test  section. 

The  remainder  of  the  gas  flow  path  consisted  of  approximately 
2.1m  of  tubing  with  the  same  diameter  as  the  test  section  followed  by 
1.4m  of  9.5mm  ID  tubing  Joined  to  the  heat  exchanger  tubing  which 
had  an  ID  of  14mm  and  was  1.9m  in  length.  A  1.3m  length  of  the  same 
diameter  tubing  continued  on  to  a  densitometer  which  had  an  inside 
diameter  of  8.3cm  and  was  61cm  long.  A  0.7m  length  of  the  14mm  ID 
tubing  continued  from  the  densitometer  to  a  plenum  chamber  of  0.14m 
ID  and  0.24m  length.  From  the  plenum  a  line  of  the  same  diameter 
as  the  entry  and  5.2m  long  ran  to  the  positive  displacement  meters 
where  the  air  was  exhausted  to  the  atmosphere.  A  flow  control  valve 
was  located  0.3  meters  after  the  plenum  chamber. 

In  low  temperature  runs,  electrical  resistance  heating  of  the  test 
section  yielded  an  axial  heat  flux  distribution  which  exponentially 
approached  a  constant  value  within  two  diameters  and  then  remained 
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constant  within  a  few  percent  until  near  the  end  of  the  section.  Tube 
wall  temperatures  were  obtained  from  premium  grade  Chromel-Alumel  thermo¬ 
couples,  0.13mm  (0.005  in.)  diameter,  which  were  spot  welded  to  the 
outside  of  the  tube.  Three  pressure  taps,  with  holes  of  about  0.30mm 
diameter,  were  used.  One  was  located  five  diameters  below  the  lower 
electrode  and  the  other  two  were  about  four  and  nine  diameters  below 
the  upper  electrode. 

The  static  pressure  at  the  test  section  inlet  was  measured  with 
a  Heise  gauge  and  the  Kulite  pressure  transducer.  Pressure  drop  was 
measured  in  steady  flow  with  an  MKS  Baratron  Model  77  differential 
pressure  meter.  Pressure  fluctuations  at  the  beginning  of  the  unheated 
entrance  section  were  measured  with  the  Kulite  transducer.  The  signals 
from  both  transducers  were  recorded  on  a  Hewlett-Packard  x-y  recorder. 
Volume  flow  rates  were  measured  with  several  Parklnson-Cowan  positive 

displacement  flow  meters  in  parallel.  The  test  section  was  completely 
enclosed  by  a  heat  shield  to  restrict  convective  air  currents  and  to 

help  stabilize  the  heat  loss  from  the  tube  to  the  environment. 

2.2.  Procedure 

To  examine  the  differences  in  heat  transfer  parameters  between 
steady  flow  and  pulsating  flow,  an  experimental  procedure  was  evolved 
to  keep  all  control  conditions  as  close  to  constant  as  possible  during 
comparison  runs. 

The  quantities  controlled  were  the  volume  flow  rate  which  deter¬ 
mines  the  mass  flow  rate  and  Reynolds  number,  the  electrical  current 
which  determines  the  heating  rate,  and  the  inlet  pressure  level.  As  a 
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consequence,  the  ges,  Reynolds  number,  Mach  number  and  non-dimensional 
heating  rate  were  held  fixed  while  only  the  pulsating  parameters, 
such  as  a,  Str,  etc.,  changed  from  zero  to  a  chosen  value.  Instru¬ 
ments  were  not  changed  during  a  set  of  runs  at  the  same  conditions, 
thereby  reducing  the  pertinent  experimental  uncertainties  to  relative 
(or  comparative)  values  rather  than  absolute  values. 

Normally  the  pulsed  run  was  conducted  first  at  the  nominal  flow 
rate,  pressure  and  heating  rate  desired.  After  the  data  were 

recorded  the  pulsed  gas  supply  was  stopped  and  the  control  settings 

were  adjusted  to  give  the  same  instrument  readings  in  steady  flow. 

Electrical  current  and  pressure  could  be  controlled  closely  but,  since 

the  flow  rate  is  deduced  from  measurements  over  a  time  interval  with 

the  positive  displacement  flowmeters,  setting  the  flow  was  an  iterative 
procedure.  The  steady  flow  rate  would  be  set  as  close  as  practical  to 

the  previous  pulsating  flow  rate,  then  after  the  data  were  recorded 
it  would  be  set  to  a  slightly  different  value  for  a  second  steady  run. 
The  latter  flow  rate  was  chosen  so  that  the  steady  flow  data  could 
be  corrected  to  the  conditions  of  the  pulsating  run  by  interpolation 
if  necessary.  Normally  the  differences  in  mass  flow  rate  were  less 
than  two  percent  so  changes  were  minimal. 

The  reproducibility  of  the  measurement  technique  was  checked  in 
two  ways.  Air  data  in  steady  flow  had  been  obtained  previously  in  two 
other  test  sections  by  Serksnis,  Taylor  and  McEllgot  [1978]  and  Pickett, 
Taylor  and  McEllgot  [1979].  It  was  found  that  each  had  a  series  of 
runs  at  Re^  near  80,000  and  various  heating  rates  so  these  were  compared 


to  present  measurements  at  the  same  conditions.  For  the  three  sets 
of  data  which  spanned  a  five  year  period,  it  was  found  that  in  the 
downstream  region  the  normalized  Nusselt  number,  Nu/(0.021Re®*®  Pr®'*), 
agreed  to  within  three  percent  at  low  heating  rates  (T^/T^  %  1.2) 
and  within  two  percent  at  higher  heating  rates  (1.4  <  T  /T,  <  1.8). 

W  D 

Secondly, the  reproducibility  of  the  present  measurements  was  tested 
at  the  end  of  the  experiments  by  duplicating  one  of  the  first  runs  with 
Be,  %  60,000,  q+  ■  0.0014  (maximum  T  /T.  %  1.5)  and  steady  conditions.  The 

1  W  D 

mass  flow  rate  could  be  reproduced  to  better  than  0.2Z,  the  test  section 
inlet  pressure  to  within  less  than  0.1Z  and  the  electrical  current 
within  the  accuracy  of  ammeter  (%  0.25Z).  The  resulting  values  of 
<Tw  -  Tfe  in)BUUt  differed  by  2.1Z,  leading  to  agreement  of  the  fully 
developed  Nusselt  numbers  within  less  than  3Z  again. 

2.3.  Range  of  variables 

Thirteen  sets  of  runs  were  conducted  with  air  (Pr  %  0.7). 

A  nominal  Reynolds  number  of  about  60,000  was  chosen  as  a  reference 

but  the  range  covered  was  19,000  <  Re  <  102,000.  Most  data  were  taken 

with  moderate  heating  (q+  %  0.0015)  with  a  few  runs  with  q+  up  to 

0.0034  (maximum  T^T^  %  2.3)  to  investigate  effects  of  property  varia- 

2 

tlon.  The  maximum  Mach  number  was  0.15  and  the  quotient  Gr^/Re^  was 
less  than  0.012  so  compressibility  and  buoyancy  effects  were  believed 
to  be  negligible.  The  operating  range  of  the  pulsating  "Gas  Booster" 
pump  was  2.1  to  3.6  Hz. 

The  response  time  of  the  test  section,  due  to  its  thermal  capacity. 


was  about  one  second  or  longer  at  the  conditions  of  the  experiments 
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Thus,  the  wall  temperature  fluctuations  were  damped  by  the  test  section 
acting  as  a  filter.  Since  the  response  time  was  not  considerably  larger 
than  the  period  of  the  forced  pulsations,  the  wall  temperature  oscillated 
slightly.  The  amplitude  was  usually  about  -  3/4*C.  with  a  few  cases 
approaching  -  1  1/2#C.  For  the  nominal  case  at  Re  3*  60,000,  the  latter 
amplitude  is  less  than  one  percent  of  the  temperature  difference 
T  -  T.  downstream.  In  data  reduction  the  mean  wall  temperature  was 

W  1 

taken  as  the  mean  of  the  two  extremes  observed  during  the  pulsations. 

For  the  pulsating  runs,  StrL  was  less  than  0.1  for  all  cases  and, 
correspondingly,  StrQ  was  0.0016  or  less.  Both  were  low  enough  to 
expect  quasi-steady  conditions  relative  to  residence  times.  The 
total  flow  length  from  the  pulsed  flow  source  to  the  exit  was  about 
20m  or  approximately  l/7th  of  the  wavelength  of  sound  at  these  fre¬ 
quencies,  so  the  llkllhood  of  significant  acoustical  resonances  is  be¬ 
lieved  to  have  been  negligible.  The  non-dimensional  frequency  or  Stokes 
parameter,  a  *  rv^2xf /v ,  ranged  from  4.1  to  7.6  and  its  thermal  counter¬ 
part  r/2irfpc/k  was  slightly  less;  if  the  flow  were  lamina^  these  mag¬ 
nitudes  would  invalidate  quasi-steady  idealizations  [Greeted  and 
Mullin,  1980],  The  measured  magnitude  of  the  pressure  pulsations, 

Ap/p  -  2(pmax  -  Pni^/CP^  +  P^).  ranged  from  0.09  to  0.29  with  the 
magnitude  of  the  bulk  velocity  pulsations  being  estimated  as  approxi¬ 
mately  half  these  values.*  Thus,  the  amplitudes  of  the  pulsations  could 
be  considered  small  to  moderate  with  none  approaching  flow  reversal. 

The  heat  transfer  data  are  tabulated  in  Appendix  B. 

_  *The  small  size  of  the  test  section  precluded  direct  measurement  of 
AV/V.  The  idea  that  AV/V  £  4  Ap/p  is  based  on  an  incompressible  flow 
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3.  RESULTS 

3. 1.  Pulsation  wave  shape.  p(8) 

The  commercial  gas  booster  pump  does  not  yield  a  pressure  fluctua¬ 
tion  that  could  be  described  as  a  pure  sine  wave,  but  typically  it  could 
be  considered  sinusoidal  as  a  first  approximation.  Figure  2a  shows 
measurements  of  the  usual  wave  form  in  these  experiments  and  Figure  2b 
shows  the  most  different  wave  form  observed.  The  latter  occurred  at 
the  slowest  frequency  used,  2.1  Hz.  Close  examination  of  the  usual 
wave  form  shows  an  almost  linear  rise  in  pressure  followed  by  an  ex¬ 
ponential  decay,  but  normally  this  decay  was  Interrupted  by  the  next 
pressure  stroke  before  it  approached  a  steady  pressure.  The  fraction 
of  time  during  which  the  pressure  was  rising  was  30  to  40  percent  in 
most  cases.  Usually,  the  decay  fraction  increased  as  the  Reynolds 
number  decreased. 


approximation  relating  Ap/p  <\J  KV2  where  K  is  an  overall  loss  coefficient. 
Prof.  Edw.  J.  Kerschen  [AME,  Univ.  Arizona,  1981]  has  suggested  an  alternate 
approach,  considering  the  system  as  a  long,  frictionless  tube  (non-steady) 
ended  with  a  restriction,  the  control  valve  (quasi-steady).  His  one¬ 
dimensional  perturbation  solution  yields  the  spatial  dependence  of  the 
velocity  pulsation  u'  in  the  pipe  as 


u*_  _  ft  sin  (k(x  -  L)/82l  +  Dk/B2cos[k(x  -  L)/S2][  g  ikM5c/6 
A  l-l  sin  [kL/S2]  +  Dk/S2  cos[kL/S2] 

where  k  «  o)/a0,  A  is  the  amplitude  of  the  imposed  pulsation,  S  »  <4  -  M2 
and  D  is  a  function  involving  the  contraction  ratio  at  the  restrictor.  He 
also  outlined  another  technique  in  terms  of  a  velocity  potential  in  order 
to  obtain  both  the  velocity  and  pressure  pulsations. 


Re  =  55, 200,  f  =  3.l5Hz,  a*6.3,  Ap/p=0.26 


One  can  apply  Che  quasi-steady  predictions  of  Mueller  [1957]  for 
rectangular  waves  to  estimate  Che  effect  of  asymmetry  in  the  wave  form. 
With  symmetrical  shapes  at  a  non-dimensional  velocity  amplitude  of 
0.2,  the  reduction  in  heat  transfer  coefficient  is  predicted  to  be  about 
one-half  percent.  With  V/V  *  0.2  for  the  rise  and  V/V  *  0.4  for  the 
decrease,  the  reduction  is  only  one  percent  so  the  sensitivity  to 
shape  or  asymmetry  is  not  expected  to  be  significant  at  moderate  amp¬ 
litudes. 

By  comparing  the  signals  measured  by  the  two  pressure  transducers, 

it  was  possible  to  estimate  the  importance  of  the  location  of  the  test 
section  transducer.  The  distance  between  the  two  transducers  was  ap¬ 
proximately  twice  the  test  section  length  and  it  included  flow  through 
an  elbow  and  a  couple  additional  fittings,  so  the  change  in  p(9)  would 
be  expected  to  be  greater  than  along  the  test  section.  For  a  given 
run  there  was  no  evident  change  in  the  shape  of  the  pressure  fluctua¬ 
tion.  The  decrease  in  the  non-dimensional  amplitude  was  approximately 
10  percent  for  Ap/p  %  0.3  and  for  Ap/p  %  0.1  it  was  up  to  40  percent. 
Therefore,  the  change  in  Ap/p  along  the  test  section  would  be  estimated 
to  be  less  than  20  percent  and  the  change  in  AV/V  would  be  expected 
to  be  less  than  10  percent.  Accordingly,  the  data  are  reported  with 
the  transducer  at  the  test  section  inlet  providing  Ap/p  for  reference. 

3.2.  Typical  thermal  entry  behavior 


section;  at  a  forcing  frequency  of  3.15  Hz,  a  was  6.3  and  Ap/p  was 
0.26,  one  of  the  higher  amplitudes  obtained.  In  the  figures,  circles 
represent  data  with  pulsating  flow. 

The  wall  temperature  measurements  are  seen  in  Figure  3  to  yield 
the  axial  profile  usually  expected  for  steady  flow  with  a  constant 
wall  heating  rate.  For  the  comparable  steady  flow  data,  the  flow  rate 
was  reproduced  to  within  1  1/2  percent  and  the  electrical  current  to 
within  0.3  percent.  The  quasi-steady  theory  predicts  only  slight 
differences  and  it  is  evident  from  even  the  raw  data  that  there  are 
not  large  differences  at  these  conditions.  In  the  thermal  entry 
most  pulsed  data  fall  between  the  two  steady  runs  as  would  be  expected 
from  the  Reynolds  numbers  if  there  were  no  effect.  Further  downstream 
the  pulsed  data  appear  slightly  lower  than  the  mean  of  the  steady 
data. 

In  order  to  eliminate  the  axial  variation  in  bulk  temperature  from 
the  comparison,  the  local  Nusselt  numbers  are  shown  in  Figure  4.  The 
differences  are  too  small  to  discern  so  it  is  concluded  that  for  these 
conditions  there  is  no  significant  effect  even  though  a  is  well  above 
two,  the  approximate  quasi-steady  limit  in  laminar  flow. 

The  slight  difference  in  Reynolds  number  may  be  treated  by  comparing 
0  8 

values  of  Nu/Re  *  .  Normalizing  then  allows  a  closer  comparison  as  in 
Figure  5.  For  this  set  of  data  the  average  difference  is  about  1/2 
to  1  percent  with  a  scatter  of  the  same  magnitude.  Most  of  the  data 
for  the  pulsed  run  are  below  those  for  the  smooth  run,  in  agreement 


Figure  3.  Comparison  of  pulsed  and  steady  measurements 
at  nominal  conditions.  Mai 7  temperatures. 


55. 2< 


Figure  4.  Comparison  of  heat  transfer  parameters  in  pulsed  and 
steady  flow  at  nominal  conditions. 


with  the  quasi-steady  predictions  of  Mueller  [1957]  and  Baird  et  al. 
[1966];  the  order  of  magnitude  of  the  reduction  is  the  same  as  predicted 
(More  exact  comparison  is  not  warranted  due  to  the  attainable  level  of 
reproducibility  mentioned  earlier.)  This  quantity, 

(nu/R«°'8)d  -  (su/Re°-8)s 

(»u/r.°-8)s 

will  be  used  in  the  rest  of  this  paper  for  comparison  of  the  pulsed 
and  steady  measurements. 

3.3.  Variation  of  parameters 

The  direct  comparisons  for  all  thirteen  sets  of  runs  yielded 

maximum  differences  for  individual  sets  ranging  from  0.7  to  7.6  percent, 

0  8 

In  most  cases  the  parameter  Nu/Re  *  was  less  for  pulsating  flow  than 
for  steady  flow.  The  effects  found  were  less  than  reported  earlier 
[Havemann  et  al.,  1956;  Romie,  1956;  Chalitbhan,  1959]  and  were  closer 
to  the  predictions  of  the  quasi-steady  analyses  (which  indicate  only 
slight  effects  at  moderate  amplitudes) . 

The  data  were  examined  for  trends  as  the  control  parameters  were 
changed.  No  consistent  variation  was  found  relative  to  non-dimensional 
amplitude;  the  largest  effects  were  found  at  low  amplitudes  rather  than 
at  high  values  as  predicted  by  quasi-steady  analyses.  Three  sets  of 
data  taken  at  higher  heating  rates,  to  see  whether  property  variation 
was  important,  showed  no  common  trend. 

The  variation  with  non-dimensional  frequency  a  is  plotted  with 
inlet  Reynolds  number  as  a  parameter  in  Figure  6.  Also  shown  on  the 
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lefC  is  the  range  o£  reproducibility  mentioned  earlier  for  the  steady 
runs;  this  range  provides  a  measure  whether  differences  should  be  con¬ 
sidered  significant.  The  symbols  identify  the  differences  calculated 
at  three  axial  locations  from  the  thermal  entry  to  downstream. 

Most  of  the  data  fall  within  the  range  of  scatter  found  in  the 

reproducibility  runs.  For  the  nominal  case  of  Re  %  6  x  10  as  well 
4  5 

as  8  x  10  and  10  ,  no  evident  trend  is  seen  as  a  varies  from  5  to  7  1/2. 
Thus,  the  present  data  for  high  Reynolds  numbers  agree  with  steady 
results  or  quasi-steady  predictions. 

The  two  points  at  Re^  %  1.9  x  10  4  do  show  an  increasing  discrepancy 
as  a  is  increased  as  do  the  few  points  at  Re^  %  4  x  104.  These  discrep¬ 
ancies  are  greater  than  the  experimental  uncertainty  in  the  direct 
comparisons.  At  a  given  non-dimensional  frequency  the  data  for  the 
lower  Reynolds  number  are  lower.  By  extrapolating  these  trends  towards 
low  a,  one  might  interpret  these  data  as  implying  that  the  limit  of 
quasi-steady  flow  may  be  given  by  a  threshold  value  of  a  which  varies 
with  the  Reynolds  number.  For  example,  from  the  figure  one  might 
estimate  a  threshold  a  of  about  2-3  for  Re^  %  1.9  x  104  and  4-5 
for  Rei%  4  x  104. 

3.4.  Discussion 

For  Re  <  4  x  104  at  a  %  5  -  7  the  present  results  showed  reductions 
in  heat  transfer  parameters  greater  than  the  one  percent  predicted 
by  the  quasi-steady  analyses  even  taking  the  steady  flow  reproducibility 


into  consideration 


The  greater  reduction  of  the  Nueselt  number  observed  at  lower  Reynolds 
numbers  and  a  %  7  can  be  considered  in  relation  to  other  criteria  sug¬ 
gested  for  turbulent  flow.  Since  convective  heat  transfer  is  dominated 
by  the  thermal  resistance  of  the  viscous  layer  (y+  30) ,  criteria 

relating  to  its  thickness  and  turbulent  bursting  phenomena  would  seem 
appropriate.  Richardson  [1967]  and  Shemer  [1981]  suggest  that  when 
the  "a.c.  boundary  layer  thickness"  becomes  comparable  to  the  "laminar" 
(linear)  sublayer  thickness  y^,  the  velocity  profiles  will  be  altered. 

Their  criterion  can  be  rewritten  as  /4irStr^/(Re^c^)  %  0.2,  an  equality 

which  is  approached  as  the  Reynolds  number  is  lowered  for  a  given 

tube  and  frequency,  but  it  is  obvious  that  if  /v/o>  »  r  (i.e.,  Str^  <<  1) 

then  /v/w  »  y  also. 

I 

Ramaprlan  and  Tu  [1980]  suggest  that  the  turbulent  transport  will 
be  modified  when  the  pulsating  frequency  is  the  same  order  as  the  tur¬ 
bulent  bursting  frequency,  fb>  Using  a  correlation  based  on  outer 
variables,  V/f^D  %  5,  by  Rao,  Narasimha  and  Badri  Narayanan  [1971], 
they  develop  the  relation  Str^  >  0.2  as  a  test  for  significant  effects. 
Chambers,  Murphy  and  McEllgot  [1982]  and  others  have  found  that  the 
bursting  frequency  scales  better  with  wall  variables  than  outer  vari¬ 
ables;  if  their  correlation  is  applied,  the  criterion  would  take  the  form 
180  StrD/(irc^ReD)  %  1.  Both  these  bursting  parameters  Increase  as 
the  Reynolds  number  decreases  for  a  given  frequency  and  tube.  Each 
of  the  above  criteria  suggest  a  modification  of  heat  transfer  param¬ 
eters  as  the  Reynolds  number  is  reduced  in  pulsating  flow.  However, 
at  Re  HS  1.9  x  10  in  the  present  test  section,  the  predicted  turbulent 


bursting  frequency  is  still  much  greater  than  the  forcing  frequency 
of  the  pulsations. 

Barnett  and  Vachon  [1970]  completed  their  analysis  with  the  con¬ 
clusion  that  the  large  changes  in  heat  transfer  noted  in  some  experi¬ 
mental  studies  must  be  due  to  a  change  in  the  basic  turbulent  exchange 
mechanisms.  Alternate  explanations  are  experimental  uncertainties 
and  systematic  error;  measurements  of  heat  transfer  in  Internal  flow 
typically  become  more  difficult  as  the  Reynolds  number  is  reduced. 

It  is  known  that  it  is  difficult  to  measure  the  flow  rate  accurately 
in  pulsating  flows  [Oppenheim  and  Chilton,  1955]  particularly  when 
a  non-linear  relation  is  involved  (this  difficulty  was  avoided  in 
the  present  study  by  the  use  of  the  positive  displacement  flow 
meters) . 

The  present  test  section  is  too  small  to  permit  direct  measure¬ 
ment  of  the  turbulence  structure  as  by  Shemer  [1981]  so  resolution 

A  4 

of  the  behavior  in  the  range  10  <  Re^  <  4  x  10  is  deferred  as  a 
topic  for  later  study.  For  the  present  we  conclude  that,  relative  to 
the  quasi-steady  predictions  in  the  low  Reynolds  number  range,  there 
is  possibly  a  small  reduction  which  increases  as  a  increases  and  as 
Re  decreases. 

A 

The  data  for  Re  >  5  x  10  essentially  confirm  the  predictions  of 
the  quasi-steady  analyses  -  that  there  is  no  significant  difference 
from  steady  flow  results  -  for  a  up  to  about  7  1/2  and  Ap/p  %  0.3  within 
the  experimental  uncertainty.  As  noted  by  Richardson  [1967],  Shemer 
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[1981]  and  others,  for  turbulent  flow  the  penetration  depth  /v/cu 

should  probably  be  evaluated  using  an  effective  viscosity  considerably 

greater  than  y  in  the  turbulent  core.  Then  the  effective  value  of 

o  ■  r  /w/v  "  would  be  reduced  compared  to  laminar  flow  at  the  same 
erf  eft 

frequency  and  tube  diameter.  This  reduction  in  a  would  be  expected  to 
become  greater  as  the  Reynolds  number  increases.  Thus,  these  data 
imply  that  a  ■  7  at  Re  >  5  %  10^  corresponds  to  ^  2  (which  would  be 

the  limit  for  quasi-steady  behavior  in  laminar  flow) . 

4.  CONCLUSIONS 

Heat  transfer  measurements  in  pulsating,  turbulent  flow  of  air 
were  compared  directly  to  results  at  the  same  values  of  the  control 
parameters  over  the  following  ranges:  4.1  <  a  <  7.6,  q+  >  0.0034 
and  19,000  Re  <  102,000.  All  pulsating  data  agreed  with  the  corres¬ 
ponding  values  for  steady  flow  within  71/2  percent. 

^  A 

For  Re  >  5  x  10  the  pulsating  measurements  essentially  confirmed 
quasi-steady  analyses,  which  predict  a  slight  reduction  in  heat  transfer 
parameters,  within  the  reproducibility  of  the  experiment.  At  lower 
Reynolds  numbers  the  reduction  was  larger  and  increased  as  a  was  in¬ 
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Appendix  A:  LITERATURE  SURVEY 
J ,  S .  Park* 

Various  studies  over  the  past  30  years  have  been  directed  at 
the  question  of  whether  pulsating  flow  can  produce  increases  in 
heat  transfer.  This  interest  is  frequently  motivated  by  the 
several  practical  situations  of  engineering  applications,  for 
example,  in  unstable  systems  of  rocket  motors,  gas  turbines,  and 
heat  exchangers  with  reciprocating  devices  or  with  compressors. 

A  survey  of  the  available  literature  for  heat  transfer  to 

pulsating  turbulent  flow  is  sumnarized  in  Table  A-l  and  Fig.  A-l. 
Generally  investigations  on  heat  transfer  to  fluids  in  pulsating 
flow  have  shown  increased  rates  of  heat  transfer  but  apparently 
conflicting  data  have  been  reported.  In  addition,  many  investi¬ 
gations  have  been  limited  in  the  range  of  variables  determining 
the  heat  transfer  phenomena  and  have  used  different  pulse  genera¬ 
tion  mechanisms,  flow  regimes,  etc. 

In  the  turbulent  flow  region,  for  example,  references 
[2,4,5,8,9,15,17,18]  show  increases  in  heat  transfer  due  to 
imparting  pulsations  to  the  fluid,  while  [1,10,11]  report  that 
these  pulsations  have  no  such  effect.  In  [3,6,7,12,13,16,19-25] 

an  Increase  in  heat  transfer  was  found  in  some  cases  and  a 
reduction  in  others,  depending  on  the  frequency,  amplitude  and 

the  Reynolds  number. 

♦Present  address.  College  of  Engineering,  Texas  A  &  M  Univ., 
College  Station,  Tex. 
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Nomenclature  for  Literature  Survey 


Symbol  Definition 


*• 

CP 

«*i 

f 

Gr 

h 

H 

k 

L 

M 

Nu 

NR 


$ 


P 

Re 

r 

Str 

t 

T 

V 


Speed  of  sound 

Specific  heat  at  constant  pressure 
Pipe  inside  diaaeter 
Pulsed  frequency,  cycles/sec 

Gras ho f  number  based  on  wall  heat  flux,  gd^/C^kT) 
Heat  transfer  coefficient 

0  Q 

Heat  transfer  parameter,  Nu/fee* 

Thermal  conductivity 
Length 
Mach  number 
Nusselt  number,  hdj/k 

Non-dimensional  parameter  for  comparisons. 


(Nu/ 


0.8, 


)„  -  <Nu/  °*8)J/(Nu/_  U*°) 


0.8, 


Re  'p  '•’"Re  's' - Re  's 

Exponent  on  heat  transfer  correlation 

Time  mean  pressure 

Peak  to  peak  amplitude  of  pulsating  pressure 
Pressure 

Non-dimensional  heat  flux,  iw/(/^CpT)in 
Reynolds  number,  /Vd^c 
Tube  radius 

Strouhal  number,  2*fdj/y 
Time 

Temperature 

Time  mean  bulk  velocity 


Nomenclature  for  Literature  Survey — Continued 


Definition _ _ 

Y  Peak  to  peak  amplitude  of  pulsating  velocity 

Y  Velocity 

Y  Non-dimensional  velocity  amplitude,  v/V 


Subscripts 
b  Bulk 

in  Inlet 

■  Mean 

p  Pulsed 

s  Steady  state 

w  Vail 

x  Axial  coordinate 

Greek  symbols 

oc  Non-dimensional  frequency,  r  \/2*f/tfr 

A  Have  length 

/ u  Dynamic  viscosity 

y  Kinematic  viscosity 

ter  Angular  frequency,  2af 


Because  the  data  reported  in  the  above  references  were 

obtained  with  different  parameters — such  as  Reynolds  number, 
frequency,  amplitude  of  pulsation,  Prandtl  number,  waveform,  etc. 

— the  apparent  conflict  of  the  results  is  not  surprising.  However, 
the  theory  of  the  effects  of  turbulent  pulsed  flow  on  heat  transfer 
is  not,  at  present,  well  known  or  understood. 

Hence,  this  paper  represents  a  contribution  toward  a 
proper  understanding  of  the  pertinent  variables  for  such  a  theory 
with  experimental  results  in  the  ranges  of  turbulent  flow  and  low 
frequency. 

In  the  present  Appendix,  a  review  on  the  effects  of  unsteady  turbu¬ 
lent  flow  generated  by  flow  pulsation  and  sound  fields  on  the  rate  of 
heat  transfer  is  presented.  The  different  studies  are  reviewed  after 
classifying  them  as  follows: 

I.  Measurements  of  heat  transfer  to  pulsating  turbulent  flow. 

A)  Water  in  steam-water  heat  exchangers  (tube) 

B)  Air  with  electrical  heating  (flat  plate) 

C)  Air  in  the  low  frequency  range  (f  ^  40Hz)  with  steam  or 
electrical  heating  (tube) 

D)  Air  in  the  high  frequency  range  (f  ^  40Hz)  with  steam  or 
electrical  heating  (tube) 

II.  Theoretical  studies  and  quasi-steady  conditions  for  pulsating, 
turbulent  flow. 

Extensive  reviews  and  bibliographies  of  the  literature  about  pulsat¬ 
ing  flows  in  general  are  also  presented  in  references  [22],  [27]  and  [28]. 
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I.  Measurements  of  Heat  Transfer  to  Pulsating,  Turbulent  Flow 

I .A.  Pulsating  Water  Flow  in  Steam-water  Heat  Exchangers 

References  [l]-[7]  have  investigated  the  effect  of 
pulsations  in  water  flow  on  the  overall  heat  transfer  coefficient 
of  a  steam-water  heat  exchanger.  Their  works  covered  the 
Reynolds  number  range  1,500-8,500  with  pulsation  frequency 

covering  the  range  0.22-16.7  Hz.  The  primary  purpose  of  these 
investigations  was  the  study  of  improving  the  overall  heat 

transfer  coefficient  by  inparting  pulsating  motion  to  the  flow  in 
the  heat  exchanger.  However,  experimental  data  pertaining  to  -this 
problem  are  clearly  Insufficient  since  they  give  conflicting  con¬ 
clusions. 

Among  those  investigations,  Karamercan  and  Gainer  [6]  and 
Herndon  et  al.  [7]  Investigated  broader  ranges  of  the  operating 
variables  than  others.  Karamercan  and  Gainer  [6]  observed  in¬ 
creases  in  the  overall  heat  transfer  coefficient  of  0-700Z  for  a 
Reynolds  number  of  1500-47,400  with  0  to  5.0  pulses/sec.  They 

j 

!  also  showed  the  highest  enhancements  in  the  heat  transfer  co¬ 

efficient  obtained  within  a  Reynolds  number  range  of  7500  to  9500. 
Herndon  et  al.  [71,  employing  0.83  to  16.7  pulses/sec  and  a 
Reynolds  number  of  6,600  to  28,000,  found  some  decreases  at  certain 
frequencies  and  Increases  at  most  frequencies  for  an  exit  pressure 
of  1  atm.  In  addition,  the  average  enhancements  observed  at 
higher  exit  pressures  (about  20-ft.  head  of  water)  were  less  than 
1.0  for  most  of  the  frequency  range  studied. 

I 
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In  general,  these  investigations  of  heat  transfer  to 
pulsatile  water  flow  in  heat  exchangers  indicate  possible  causes 
of  enhancement  as  follows: 

1)  to  increase  the  level  of  turbulence  in  the  pulsed  stream 

2)  to  cavitate  the  fluid  next  to  the  tube  wall 

3)  to  produce  periodic  reversals  in  the  pressure  gradient. 
Karaaercan  and  Gainer  [6]  state  that  periodic  reversals  in  the 
direction  of  flow  seemed  to  be  the  important  parameter  in  the 
enhancement  of  heat  transfer  in  pulsating  flow  because  flow 
reversal  increases  cavitation  and  also  promote  a  higher  level  of 
turbulence  in  the  fluid.  However,  the  manner  in  which  these 
mechanisms  Improve  the  convective  heat  transfer  coefficient  is  still 
somewhat  vague. 

I.B.  Pulsating  Air  Flow  Over  a  Flat  Plate 

References  [8]-[ll]  are  concerned  with  turbulent  heat 

transfer  from  flat  plates  immersed  in  oscillating  flows  with 
frequencies  from  3.0  to  680  Hz  and  Reynolds  numbers  from  10,000  to 
48,000  as  shown  in  Table  A-l.  Investigations  by  [8]  and  [9]  resulted  in 
increases  in  Nusselt  number  of  30  or  65Z.  In  contrast  Failer  [10] 
and  Miller  [11]  found  practically  no  difference  in  heat  transfer 
rates  obtained  with  and  without  pulsations  in  the  frequency 
range  3.0-200  cps.  Unfortunately,  those  investigations  are 
limited  to  a  range  of  Reynolds  numbers  from  10,000  to  48,000, 
approaching  the  transition  region  in  such  flows. 


Table  A-l.  Summary  of  data  on  heat  transfer  to  pulsating  turbulent  flow. 


(continued) 


Table  A-l  (continued) 


I.C.  Pulsating  Air  Flow  in  Tube  in  the  Frequency  Range  0-40  Hz 
The  results  of  measurements  of  heat  transfer  to  pulsating 
air  flow  through  a  tube  in  the  low  frequency  range  are  shown 
in  Fig.  A-2  and  Fig.  A-3.  Half  [12]  [14]  [15]  of  those  experiments 
used  steam-to-air  heat  exchangers  which  do  not  normally  permit 
precise  comparison,  and  the  other  experiments  [131,  [16],  [25] 
used  the  tubes  with  electrical  heating. 

By  using  a  tube  heated  with  steam,  Haveman  et  al.  [12] 

(f-5-40  Hz,  Re*6000-25,000)  reported  changes  of  heat  transfer 
parameters  from  -40%  to  +40%;  Chalitbhan  [15]  (f”0-15Hz,  Re- 7 000- 
200,000)  always  showed  an  increase  of  heat  transfer,  as  high  as 
100%  at  Re  * 10, 000-50, 000  and  around  20%  at  Re*160, 000-200, 000; 
and  Mueller  [14]  (f-0. 038-0. 248  Hz,  Re-53,000-76,000)  showed  the 
average  Nusselt  number  to  be  about  0-20%  less  than  the  corres¬ 
ponding  steady  flow,  theoretically  and  experimentally.  Because 
no  information  was  given  on  the  amplitude  of  the  fluctuations 
and  substantial  experimental  uncertainties  were  involved  (generally, 
more  than  ±10-30%  for  steady  flow) ,  it  is  very  difficult  to 
compare  these  results  with  other  references. 

Particularly,  Chalitbhan  calculated  Nusselt  numbers  by  modi¬ 
fying  the  Martlnelll  equation  [42],  and  computed  the  Nusselt 
number  of  the  pulsating  flow  by  using  the  ratio  of  the  pressure 
drop  for  pulsating  flow  to  that  for  non-pulsating  flow.  He  used 
a  water-filled  manometer  for  the  pressure  drop  of  the  test 
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Fig.  A-2.  The  effect  of  air  velocity  on  heat  transfer  to  pulsating,  turbulent  flow 
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Frequency,  Hi 

Fig.  A-3.  The  effect  of  the  frequency  of  induced  fluctuations  on  heat 
transfer  to  pulsating,  turbulent  air  flow. 


section.  So  his  results  nay  have  a  large  error  as  the  values  of 
pressure  drop  for  pulsating  flow  cannot  be  measured  accurately  by 
inferential  meter  readies  [43]  [44]. 

With  electrical  heating,  Romie  [13],  Mamayev  et  al.  [16]  and 
Creff  et  al.  [25]  showed  that  heat  transfer  parameters  of  the 
pulsed  flow  are  close  (±202)  to  those  obtained  in  a  steady  flow, 
except  for  the  resonance  results  of  Ref.  [25]  (Fig,  3)* 

Romie  [13]  measured  heat  transfer  data  at  a  position  12 
diameters  downstream  from  the  heated  section  entrance  (Re-5000, 

A 

v 

f *3. 3-133  Hz,  *0. 5-2.0)  and  showed  the  ratio  of  the  Nusselt 
numbers  of  the  pulsed  flow  to  those  of  unpulsed  flow  increases 

A  . 

above  unity  with  increasing  frequency  (  V/V>  1.0)  then  decreases  to 
a  value  below  unity,  and  then,  at  still  higher  frequencies, 
increases  above  unity.  But  as  explained  by  Bogdanoff  [22]  for 
the  high  frequencies  (f-75,  133  Hz),  the  oscillation  amplitudes 
measured  by  a  hot-wire  would  not  be  representative  of  conditions  at 
the  point  where  the  heat  transfer  measurement  was  taken,  because 
becomes  much  shorter  than  the  duct  length  and,  hence,  the  velocity 
amplitudes  at  a  point  may  not  be  representative  of  conditions 
existing  all  along  the  heated  section. 

Mamayev,  et  al.  [16]  carried  out  experiments  (f-0.5-25  Hz, 
Re-540-11,000)  and  showed  that  the  heat  transfer  coefficient 
averaged  over  the  pipe  length,  is  higher  (0-152)  than  in  an 
equivalent  non-pulsed  flow  except  for  f«20  Hz  as  shown  in  Figs.  A- 2 
and  A-3.  The  amplitude  of  fluctuations  was  not  recorded.  Their 


effects  of  air  velocity  and  the  frequency  on  heat  transfer  to 
pulsed  flow  at  Re«4, 000-8, 800  can  be  summarized  as  follows: 

I.  The  effect  of  air  velocity  on  heat  transfer  to  pulsed  flow 

A)  The  critical  values  of  Reynolds  number,  reflecting 
changes  in  flow  modes,  depend  on  f  and  shift  toward  lower 

values  with  an  increase  in  f.  This  effect  was  not  observed 
at  f£2  Hz. 

B)  At  Re>  the  ratio  of  heat  transfer  coefficients  drops 

exponentially  and  steeply  with  Re  for  f-5,  12,  20  Hz. 

C)  The  ratio  of  heat  transfer  coefficients  tends  to  unity  at 
f-0.5  to  12  Hz,  while  at  f^  20  Hz  the  pulsations  have  a 
negative  effect  on  the  heat  transfer  rate  in  turbulent 
flow. 

II.  The  effect  of  the  frequency  on  heat  transfer 

A)  At  Re-3000  to  4000,  the  ratio  of  heat  transfer  coefficients 
Increases  above  f-5  Hz,  and  then  asymptotically  tends  to 
unity. 

B)  At  Re-5000-8800,  this  ratio  has  a  peak  at  f-0.5  to  1.0  Hz 
and  then  decreases  to  unity. 

C)  No  effect  of  flow  fluctuations  is  observed  at  f«15  Hz, 
while  at  f-15  to  24  Hz  the  heat  transfer  coefficients  are 
smaller  than  for  steady  flow. 

As  shown  in  Fig.  A-3,  Romie  and  Mamayeu  et  al.  showed  no  effects 
of  flow  fluctuation  at  a  frequency  near  IS  Hz.  Creff  et  al.  [25], 
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(mean  Re^lO^  -2x10^,  f»0-36  Hz,  —  ■  0.36),  studied  the  frequency 

V 

influence  on  the  local  heat  transfer  rates,  and  especially  those 
due  to  acoustic  resonance  frequencies  of  the  pipe.  They  showed 
that  the  local  Nusselt  numbers  are  close  (+5.0~-20Z)  to  those 
of  the  unpulsed  flow  at  the  same  mean  flow  rate  except  near  the 
resonance.  For  the  resonance  modes  heat  transfer  increases  about 
75%  at  antinodes  of  the  amplitude  of  the  gas  velocity  as  shown  in 
Fig.  A-3.  Heat  loss  is  less  than  or  equal  to  10%  of  total  heat  flux. 

l.D.  Pulsating  Air  Flow  in  Tube  in  the  High  Frequency  Range,  f  >  40  Hz 

Experimental  Investigations  in  the  high  frequency  ranges 
which  Involved  resonances  [26]  of  the  tube  system  were  made  by 
using  a  steam  heated  tube  with  acoustic  vibrations  (siren  or  horn) 
[17-20,  22]  and  by  using  an  electrically  heated  tube  with 
oscillations  generated  by  a  rotating  valve  [21,  23,  24]. 

Most  of  these  results,  generally,  agreed  that  the  effect  of 
resonant  fluctuations  on  heat  transfer  was  found  much  greater  than 
that  of  non-resonant  pulsations.  The  effect  increases  as  the 
amplitude  of  pulsations  increases. 

However,  the  effects  of  resonant  fluctuations  in  heat  transfer 
are  uncertain  because  an  increase  in  heat  transfer  was  observed  in 
some  cases  and  a  reduction  in  others,  depending  on  the  Reynolds 
number,  node  or  antinode  of  standing-wave  velocity  along  the 
length  of  pipe,  pulsation  amplitude,  frequency,  and  so  on. 
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Hwu  and  Lemlich  [17,  18]  shoved  an  increase  in  the  heat  trans¬ 
fer  coefficient  up  to  27%  in  the  turbulent  region,  at  frequencies 
which  were  resonances  of  the  tube  system  allowing  large  amplitude 

oscillations  to  occur  (Re*560-5,900,  f«198,  256,  322  Hz) .  For 
pressure  amplitude  measurement  a  water  manometer  was  used,  so  the 

amplitude  data  may  not  be  accurate,  as  mentioned  in  Ref.  [43,  44]. 

Jackson  et  al.  [19,  20]  studied  air  flow  in  a  steam  heated 
tube  (Re-2000-200,000,  f*90-356  Hz)  at  resonant  frequencies.  To 
measure  the  intensity  of  the  oscillations  induced  by  a  horn,  a 
microphone  mounted  on  a  rod  was  inserted  into  the  duct  from  the 
upstream  end.  Typical  heat  transfer  data  are  reproduced  in 
Figs.  A-4  -  A-7 ,  and  are  discussed  by  Bogdanoff  [22]  for  a  frequency 
of  about  220  cps  as  follows: 

1)  At  Reynolds  numbers  of  43,000  and  above,  the 
data  behaves  as  exemplified  by  Fig.  A-6,  i.e., 
the  general  effect  of  oscillations  is  to  reduce 
the  heat  transfer  rates,  the  largest  reductions 
appearing  near  the  velocity  antinodes,  and 
relatively  little  effect  in  the  regions  of  the 
velocity  node. 

2)  As  Reynolds  number  is  decreased  into  the  range 
22,800-33,000  the  effects  of  heat  transfer 
become  small  and  irregular  although  the  sound 
pressure  levels  are  diminished. 

3)  On  further  decrease  of  the  Reynolds  number  into 


A-19 

the  range  11,600-16,000,  there  ia  again,  a  strong 
effect  of  the  oscillation  on  heat  transfer  as  shown 

in  Fig.  A-5,  i.e.,  increases  in  heat  transfer  at 

the  velocity  antinodes,  decreases  at  the  velocity  nodes 

and  an  overall  increase  in  heat  transfer. 

Bogdanoff  [22]  investigated  air  flow  in  a  steam  heated  pipe  at 
the  resonant  frequency  which  was  induced  by  a  siren  wheel. 

Measurement  of  the  amplitude  of  the  pressure  fluctuations  was 
taken  with  a  pressure  transducer  just  upstream  of  the  siren  wheel 
which  was  at  the  end  of  tube.  The  experimental  results  for  7  runs 
at  the  9th  harmonic  and  1  run  at  the  13th  harmonic  are  reproduced 
in  Figs.  A-8  -  A-15. 

The  quantity  uA/jj  .  the  Instantaneous  velocity  which  would  be  required 

*•  mm 

to  cause  a  pressure  fluctuation  p/p  ,  was  computed  from  a  simple-wave 
formula  for  acoustical  waves  in  a  duct  without  friction  losses 
[26]  , 

UA/jj  -  pN/_  *  1/rM  (A.l) 

where  the  pertinent  nomenclature  may  be  summarized  as: 

S*  *  wD /|Ja,  Strouhal  number  upstream  of  the  heated  section 
PN  ■  root  mean  square  pressure  fluctuation 

A# 

U  “  r.n.s.  velocity  fluctuation  at  pipe  center-line 

A 

0  •  peak-to-peak  velocity  at  pipe  center-line,  or  in  free 
stream 

§,  U  ■  mean  pressure,  velocity. 

A  mm 

From  the  assumption  of  sinusoidal  pressure  waves,  U^/U  is  roughly 


0*/jj=2  VT?a/5=2.83  5x/j 


(A. 2) 
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Ma  «  .0309 

Gat- Air 

So  «  6.34 

A  • 34.7  ptia 

P„/P  •  0199 

Pipe  LO.  *  1.3  in 

9  th  Harmonic 

Uj/U  *  .433 

Rea  •  100,200 

Steady  flow 


Distance  from  entrance  of  heated  section,  dias. 


I  234  5  6  7  8  9  to  it  12  13  14  15  16  17  18  19  20  21  22  23  24  23 

Chamber  no. 

Figure  A-8.  Effect  of  pulsations  on  local  Nusselt  numbers  for  p/-  *  0.0196 
from  Bogdanoff  122]. 
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Distance  from  entrance  of  heated  section,  dias. 


1234  5  6  7  8  9  iO  II'  12  13  14  15  16  17  18  19  20  21  22  23  24  25 

Chamber  no. 

Figure  A-9.  Effect  of  pulsations  on  local  Nusselt  numbers  for  p/-  »  0.0343 
from  Bogdanoff  [22]. 


His  results  can  be  summarized  as  follows: 

1)  For  all  runs,  Re  %  105,  and  M  -  0.03,  S  %  6.4  for  all 

a 

9th  harmonic  runs  and  Sa-9.18  for  the  13th  harmonic  run. 

From  Fig.  A-8  we  see  that  oscillations  have  little  effect  on  heat 
transfer  at  UA/jj —0.4. 

2)  At  U^/^%0.8,  the  oscillations  produce  a  noticeable  effect 
on  the  heat  transfer.  The  maxima  of  heat  transfer  are  downstream 
of  the  velocity  nodes,  the  minima  downstream  of  antinodes  and  there 
is  little  overall  change  of  heat  transfer. 

3)  As  the  oscillating  amplitude  increases  (  0A/fj  a*  1.2,  1.4, 
1.56),  the  amplitude  of  the  variation  of  heat  transfer  increases, 
and  the  maxima  and  minima  move  downstream,  but  the  overall  heat 
transfer  changes  are  small. 

4)  Increasing  the  oscillation  amplitude  still  further 

(  0A/B  -  1.7,  2  .3)  moves  the  maxima  and  minima  even  further 
downstream  (the  maxima  is  near  the  velocity  antinode,  and  the 
minima  are  slightly  downstream  of  velocity  nodes),  increases  the 
amplitude  of  the  heat  transfer  fluctuations,  and  produces  sub¬ 
stantial  overall  increases  in  heat  transfer  rates  (especially  for 

Vo  *  2*3>* 

5)  The  test  taken  at  the  13th  harmonic  shows  fluctuations  of 
heat  transfer  whose  maxima  and  minima  are  located  similarly  to 
those  at  the  9th  harmonic  for  comparable  UA/^j  values 

(  UA/o  *  1.9)  but  shows  considerably  smaller  overall  Increases  in 


heat  transfer 


Koshkin  et  al.  [21]  and  Galitseyskiy  et  al.  [23,  24]  investi¬ 
gated  heat  transfer  to  pulsating  air  flow  in  an  electrically 
heated  tube.  Local  heat  transfer  measurements  were  taken  at  various 
points  over  a  heated  section  (Re  *  10^-10^,  p  -  72-290  psia, 

A 

P 

Tw/Tb  «  1. 2-1.6  [23,  24]  or  1.0  [21],  —  -  0-0.25). 

P 

Oscillations  were  generated  by  a  rotating  valve  upstream  of  the 
heated  section.  The  relative  amplitude  of  pressure  pulsations  was 
taken  at  the  experimental  tube  inlet.  Typical  heat  transfer  data 
are  reproduced  in  Fig.  A— 16.  These  results  can  be  summarized  as 
follows: 

(1)  Resonance  pressure  pulsations  of  the  fluid  in  a  pipe 
appreciably  affect  the  heat  transfer  near  standing-wave 
velocity  maxima;  the  effect  increases  with  an  increase  in 
the  amplitude  of  the  pressure  pulsation. 

(2)  In  the  experiments  the  heat  transfer  coefficient  in  the 
first  standing-wave  velocity  maxima  was  2-3  times  greater 
than  that  for  steady  state. 

(3)  As  a  result  of  dissipation  of  the  pulsation  energy,  the 
heat  transfer  coefficient  decreases  along  the  pipe,  i.e., 
the  closer  the  velocity  maximum  to  the  pipe  inlet  the 
higher  the  heat  transfer. 

(4)  The  distribution  of  the  local  heat  transfer  coefficient 
along  the  length  of  the  pipe  is  similar  to  the  kinetic 


Figure  A-16.  Heat  transfer  distribution  (K  *  Nup/Nus)  along  the 
experimental  tube  with  n  ■  3  for  various  values  of 
(7V  K»  -  'O4  T  105;  a  -  (*£),  ■  0.225;  b  -  (4)„  -  0.184; 

c-f—')  *  0.11,  d  -  (~z\  *  0.090,  e-f^\  =  0.054. 

\p/o  \  p  a>  \  p  a> 

(from  Koshkin  et  al.  [21]} 
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energy  distribucion  along  che  standing  wave. 

(5)  The  number  of  heat  transfer  maxima  and  minima  is  determined 
by  the  number  of  the  resonant  harmonics. 

(6)  The  effects  of  the  Reynolds  number  in  the  range  of  10^-10^ 
and  of  the  number  of  the  resonance  harmonic  on  the 
relative  heat  transfer  in  the  experiments  are  insignificant 
as  under  steady-state  conditions,  i.e.,  lie  within  the 
limits  of  experimental  accuracy  (+  10%) . 

II.  Theoretical  Studies  and  Quasi-steady  Conditions  for  Pulsating 
Turbulent  Flow 

Although  the  effects  of  the  flow  pulsation  on  heat  transfer 
characteristics  for  turbulent  flow  have  been  studied  over  the  past 
40  years,  there  are  few  theoretical  analyses  in  the  literature 
[14,29,31]  for  turbulent  pulsating  flow.  The  present  work  is 
limited  to  demonstrate  simply  the  results  of  the  theoretical 
analyses  rather  than  a  comprehensive  explanation  of  these  studies. 

Barnett  and  Vachon  [29]  presented  an  analysis  for  the  turbulent 
fully  developed  flow  of  a  fluid  in  a  tube  undergoing  harmonic 
oscillations  parallel  to  its  centerline  with  the  governing  param¬ 
eters:  Reynolds  and  Prandtl  numbers,  non-dimensional  frequency, 
and  the  dimensionless  amplitude  of  vibration.  By  assuming  the 
turbulent  diffusion  of  momentum  and  energy  were  unaffected  by 
this  motion,  they  predicted  significant  Increases  in  heat  transfer 
coefficients  only  at  low  frequencies  and  for  large  amplitudes. 
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The  Nusselt  number  decreases  at  high  frequencies  with  respect  to 
the  corresponding  steady  flow  value.  The  effects  of  pulsating 
flow  on  heat  transfer  are  amplified  with  Prandtl  numbers  below 
unity. 

References  [3,  4,  5,  14,  22,  30,  31,  32]  discussed  the  effects 
of  pulsatile  flow  to  heat  transfer,  and  some  provided  results  for 

quasi-steady  conditions;  i.e.,  frequency  low  enough  that  heat 
transfer  coefficient  at  any  instant  in  pulsating  flow  can  be  pre¬ 
dicted  by  the  usual  steady  state  correlations. 

Park,  Taylor  and  McEligot  [33]  summarized  the  comments  of 
Smolder en  [34]  about  the  governing  parameters  of  non-steady  flow 
as  follows: 


1)  Strouhal  number:  StrL  ■ 


(A.3) 


2)  For  acoustic  disturbances,  the  ratio  of  the  length 
of  the  flow  path  to  the  wave  length  of  sound  waves: 


L/a  •  —  ■  M  .  StrT 
ao  L 


(A. 4) 


3)  For  consideration  of  viscous  effects,  non- 
dimensional  frequency: 


oc  -  L  V 2  7t  f  /y  (-'jReL  .  StrL  ) 


(A.  5) 


If  StrL  4C1,  the  unsteady  term  in  the  governing  equation  becomes 
negligible  relative  to  convective  terms,  and  the  flow  can  be 
treated  as  quasi-steady.  If  ^  and  at  are  small,  the  flow  is 
considered  as  quasi-steady  for  acoustical  and  viscous  effects. 


respectively. 
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Mueller  [14]  showed  analytical  and  experimental  results  of 
pulsating  turbulent  pipe-flow  heat  transfer  in  a  quasi-steady 
condition,  for  quasi-steady  flows  which  are  non-reversing,  his 
analysis  reported  a  slightly  lower  average  Nusselt  number  must 
be  expected  than  for  steady  flow.  Experimental  results 
(Re*5 3 ,000-76,000,  f “0.038-0. 248  Hz,  steam  heated  tube)  showed 
the  average  Nusselt  number  to  be  less,  0-20?.,  than  the  corresponding 
steady-flow  Nusselt  number. 

Lemlich  [3]  showed  that  for  fully  developed  flow  with  constant 
properties,  if  h8  Vn,  the  improvement  ratio  or  enhancement  due 
to  pulsation  would  be 

_  (2X)11*1!  v®(*ut)  d(»t) 

hn/v,  *  —  -  (A. 6) 

^  £l  ,(“Jt)  d<"t)]° 

By  similar  method,  Baird  et  al.  [4,  5]  and  Bogdanoff  [22] 
showed  that  the  tlme-aveTage  ratio  of  heat  transfer  can  be  shown 


to  be. 


5“  1  f27l|  V  l°»8 

P/Nug  *  SJ  I1  *  s4iwt  I  d(«t) 


(A.  7) 


0  8 

with  the  assumptions  of  hs  /w  V  *  for  steady  turbulent  flow  in 
a  tube  and  a  sinusoidal  variation  of  velocity. 


V  -  V  (1  +  -y  sin  WJC) 


(A. 8) 


Typical  results  of  Eq.  A. 7  which  are  Independent  of  frequency 

are  shown  in  Fig.  A-17.  Experimental  results  by  Baird  et  al.  [4,  5] 
showed  that  h^/hs  does  not  fall  significantly  below  unity  at  low 
values  of  V/V  ,  and  that  the  quasi-steady  flow  theory,  while  not 


<>  <> 


/v  <2.0 

/V  >2.0 


t  no  reverse  flow 
t  reversed  flow  occurs 
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completely  satisfactory,  does  predict  approximately  the  effects  of 
flow  pulsation  at  high  values  of  V/V. 

It  is  shown  in  Fig.  A-17  that  significant  improvements  in  heat 
transfer  are  not  obtained  until  V/^f  exceeds  approximately  3.0. 

Many  papers  [4,  14,  28,  31,  32,  35,  36]  discussed  different 
conditions  for  the  quasi-steady  state  approximation  to  be  valid, 

e.g.,  oC  <7.4  [4]  ,  oCZ^0.1  Bejj  [35,36]^ 

III.  Summary  of  Literature  Survey 

The  results  presented  in  this  literature  survey  show  that 
different  studies  of  heat  transfer  in  pulsed  flow  have  provided 
their  contributions  by  showing  the  existence  of  specific  phenomena 
appropriate  to  those  flows.  Those  contributions  have  shown  that 
the  heat  transfer  rate  increases  or  decreases  in  comparison  to 
non-pulsed  flow  for  various  experiment  conditions,  by  variation 
of  the  following  quantities: 

1)  frequency 

2)  amplitude  of  oscillation 

3)  mean  Reynolds  number 

4)  mode  of  generation  of  pulsations 

5)  acoustic  resonance  of  the  system 

6)  antinode  or  node  of  the  flow  velocity 

7)  and  so  on. 

It  is  not  immediately  apparent  which  of  the  parameters  are 
most  important  in  determining  the  nature  of  the  effects  of  the 


pulsations  on  heat  transfer.  It  is  likely  that  different  para¬ 
meters  will  dominate  in  different  ranges.  The  analyses  generally 
predict  only  a  slight  modification  of  heat  transfer  parameters  in 
pulsating  turbulent  flow  whereas  experiments  have  found  larger 
effects. 

Therefore,  it  is  important  to  measure  the  heat  transfer  param¬ 
eters  in  typical  flows  where  data  are  not  available  in  order  to 
test  the  analytical  predictions  for  normal  operating  conditions. 
Also,  we  should  look  for  the  parameters  which  most  influence  the 
heat  transfer  in  pulsating  turbulent  flows.  This  is  the  task 
recommended  for  the  current  and  future  studies. 
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APPENDIX  B:  EXPERIMENTAL  RESULTS 


Table  B-l.  Summary  of  experimental  conditions 


Run  Frequency  (Hz)  a 

Ap/p 

Ret  x  10~4 

<vv 

805 

2.10 

5.35 

0.13 

6.14 

1.53 

806 

— 

— 

— 

6.00 

1.54 

807 

— 

— 

— 

6.18 

1.52 

808 

3.15 

6.18 

0.26 

5.52 

1.51 

809 

— 

— 

— 

5.43 

1.52 

810 

— 

— 

— 

5.59 

1.51 

811 

3.56 

4.31 

0.20 

1.88 

1.50 

812 

— 

— 

— 

1.91 

1.49 

813 

— 

— 

— 

1.90 

1.49 

814 

2.84 

7.08 

0.09 

1.93 

1.51 

815 

— 

— 

— 

1.94 

1.48 

816 

— 

— 

— 

1.94 

1.48 

817 

2.80 

7.17 

0.15 

3.76 

1.51 

818 

— 

— 

— 

3.84 

1.49 

819 

— 

•m  -i 

— 

3.79 

1.49 

820 

2.71 

7.18 

0.26 

5.78 

1.88 

821 

— 

— 

— 

5.82 

1.91 

822 

— 

— 

— 

5.90 

1.89 

823 

2.56 

7.03 

0.29 

9.65 

1.49 

824 

— 

— 

— 

10.22 

1.48 

825 

— 

— 

■— 

9.93 

1.49 

826 

2.81 

6.65 

0.35 

7.69 

1.49 

827 

— 

— 

— 

7.85 

1.49 

828 

— 

— 

7.64 

1.49 

829 

2.54 

7.53 

0.28 

9.50 

2.20 

830 

— 

— 

— 

9.61 

2.18 

831 

— 

— 

— 

9.43 

2.19 

832 

2.84 

7.41 

0.27 

7.55 

2.28 

833 

— 

— 

— 

7.48 

2.27 

834 

— 

— 

— 

7.60 

2.26 

835 

2.94 

6.68 

0.28 

7.74 

1.27 

836 

— 

— 

— 

7.82 

1.26 

837 

— 

— 

— 

7.72 

1.26 
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Table  B-2.  Tabulated  data 


The  headings  used  In  the  following  listings  of  the  heated  flow  data 
and  their  definitions  are  below. 


PR,  IN 
GR/RESQ 


MACH (2) 

MACH (16) 

T.SURR 

TC 

X/D 

TW 

TW/TB 

HL/QGAS 

QGAS 


PT 

TB 

PRESS  DEFECT 


Definition 
Inlet  mixer  temperature 
Calculated  exit  temperature 
Alternating  current 
Voltage  drop  between  voltage  taps 
Inlet  Prandtl  number 

Ratio  of  Grashof  number  to  the  square  of  the 
Reynolds  number 

Mach  number  at  thermocouple  2 

Mach  number  at  thermocouple  16 

Temperature  of  surroundings 

Thermocouple  number 

Axial  position,  corresponds  to  x/D  in  text 
Inside  tube  wall  temperature,  °F 
Wall-to-bulk  temperature  ratio 
Ratio  of  heat  loss  to  heat  flux  to  gas 
Heat  flux  to  gas 

Non-dimensional  turbulent  heat  flux  parameter. 
Corresponds  to  q+  in  text 

Pressure  tap:  1-near  inlet,  2-near  exit 

Bulk  static  temperature 
2 

°i8c(pi  "  p)/G 
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RUN  824H*  DATE  8/06/61*  GAS  AIR ( STEADY)  >  MOLECULAR  WT.  >  28.97 

TIN  ■  79.5  F,  TOUT  ■  219.7  F#  MASS  FLOW  RATE  ■  68.9  LB/HR*  I  -  106.0  AMPS*  E  •  6.966  VOLTS 
PR* IN  •  .719*  GR/RESO  -  .226E-02*  MACH (2 )  -  .129*  MACH<16)  «  .142*  T*SURR  >  112.0  F 
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RUN  826H»  DATE  8/06/81  *  GAS  AIR(PULSED)  #  MOLECULAR  UT.  >  28.97 

TIN  •  88.8  F •  TOUT  -  251.8  F#  MASS  FLOW  RATE  -  52.5  LB/ HR#  I  -  100.0  AMPS#  E  -  6.550  VOLTS 
PR#  IN  >  .717#  GR/RESQ  >  .191E-02#  MACH<2)  •  .126#  MACHU6J  •  .143#  T#SURR  >  108.5  F 
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RUN  827H*  DATE  8/08/81  *  GAS  AIR ( STEADY )  *  MOLECULAR  UT.  ■  28.97 

TIN  -  86.6  F>  TOUT  >  246.9  F*  MASS  FLOW  RATE  »  53.4  LB/HR#  I  •  100.0  AMPS*  E  >  6.515  VOLTS 
PR*IN  -  .716*  GR/RESQ  >  .205E-02*  MACH( 2 )  *  .125*  MACH(16>  •  .140*  T*SURR  >  106.7  F 
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RUN  8 30H*  DATE  8/09/81  *  GAS  AIR ( STEADY)  *  MOLECULAR  WT.  •  28.97 

TIN  -  88.8  F#  TOUT  -  458.7  ft  MASS  FLOW  RATE  >  65.7  L  B  /  HR  >  1  -  166.8  AMPS*  E  >  11.250  VOLTS 
PR»IN  •  .717*  GR/RESQ  >  .721E-02,  MACHC2)  •  .112*  MACH(16)  -  .140*  T*SURR  •  196.0  F 
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RUN  8 38H#  DATE  8/10/81  #  GAS  AIR(PULSED)  #  MOLECULAR  WT.  >  28.97 

TIN  •  82.1  F»  TOUT  -  162.9  F#  MASS  FLOW  RATE  •  24.3  LB/HR*  1  •  48.6  AMPS#  E  »  3.232  VOLTS 
PR# IN  >  .718#  GR/RESO  •  .259E-02#  MACH ( 2 )  •  .054#  MACHC16)  -  .056#  T#SURR  •  88.3  F 
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